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ABSTRACT 

We report the ATCA and ROSAT detection of Supernova Remnant (SNR) J0529-6653 
in the Large Magellanic Cloud (LMC) which is positioned in the projected vicinity 
of the known radio pulsar PSRB0529-66. In the radio-continuum frequencies, this 
LMC object follows a typical SNR structure of a shell morphology with brightened 
regions in the south-west. It exhibits an almost circular shape of D=33x31 pc (1 pc 
uncertainty in each direction) and radio spectral index of a=-0.68±0.03 - typical 
for mid-age SNRs. We also report detection of polarised regions with a peak value 
of ~17%±7% at 6 cm. An investigation of ROSAT images produced from merged 
PSPC data reveals the presence of extended X-ray emission coincident with the radio 
emission of the SNR. In X-rays, the brightest part is in the north-east. We discuss 
various scenarios in regards to the SNR-PSR association with emphasis on the large 
age difference, lack of a pulsar trail and no prominent point-like radio or X-ray source. 

Key words: supernova remnants - pulsars - Large Magellanic Cloud - 
SNRJ0529-6653. 



1 INTRODUCTION 

The Large Magellanic Cloud (LMC) provides an excellent 
laboratory to study supernova remnants (SNRs) at a known 
distance of 50 kpc ( di Benedetto||2008 ). The line of sight to 
the LMC lies well away from the Galactic plane, minimising 
the obscuration and confusion from the foreground gas, dust 
and stars. 

A distinguishing characteristic of SNRs in radio wave- 
lengths is their predominantly non-thermal continuum emis- 
sion. Generally, SNRs display a radio spectral index of 
a ^ —0.5 (defined hy S oiv^)^ although a may vary signifi- 
cantly, as there exists a wide variety of types of SNRs in dif- 
ferent environments and stages of evolution (Filipovic et al. 



1998). For example, younger remnants can have a spectral 
index of a ^ —0.8, while older SNRs and Pulsar Wind Nebu- 
lae (PWN) tend to have flatter radio spectra with a ^ —0.2. 
SNRs have a great impact on the physical properties, struc- 



ture, and evolution of the interstellar medium (ISM). Con- 
versely, the interstellar environments in which SNRs reside 
will heavily affect the remnants' evolution. 

Type II supernovae result from the core collapse of 
massive stars with initial masses greater than ^ 8 ± IMq 



(Smartt 2009), and may leave behind compact central ob- 



jects, such as neutron stars that may be observable as pul- 
sars. However, not many SNRs from Type II supernovae are 
observed to host pulsars. Among the 51 confirmed and 25 
possible candidat e SNRs in the LMC ( [Klimek et al.||2010 



Desai et "aL]|2010 Filipovic et al. in prep), only three (N49, 



30 Dor B, B0540-693) have documented associations with 
a pulsar fManchester et al. 2006| . There are over 70 well 
studied SNRs in the Magellanic Clouds (MCs) as well as 
19 detected pulsars. Ridley & Lorimer (2010) modelled the 



potentially observable radio pulsars in the MCs and pre- 
dicted som e 1.79x10^ p ulsars. The Milky Way (MW) has 
274 SNRs ( |Green|[2009t and -1900 pulsars. Therefore, the 
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observed SNR/PSR ratio in the MCs is much less than in 
the MW (~l/20). However, if the pulsar surveys are sensi- 
tivity limited, then the expected ratio of the mean distances 
^(10/50)^ would result in the two fractions being similar. 
This rarity of pulsar-SNR connections may be attributed 
to the fact that many young neutron stars in SNRs exhibit 



different properties from the traditional radio pulsars ( Got- 
[thelf Vasisht'2QQ0 ). Radio pulsars remain detectable well 
after their SNRs have merged into the ISM; therefore, most 
pulsars are not in SNRs and their general properties may be 
different from those young ones in SNRs. 

PSRB0529-66 in the LMC was the first extra-galactic 
pulsar discovered. It was found by [McCulloch et al.| 
(1983) using the 64-m radio telescope at the Parkes 
Observatory. The pulsar was studied by [Costa et al.| 
( 1991 ) at a frequency of 600 MHz with the intent of 
using a polarimeter to collect integrated profiles, rota- 



tion and dispersion measures. Crawford et al.| ( |200l| ) used 
the Parkes telescope at 20 cm and improved the radio 
positional accuracy to RA (J2000)=05^2 9^50.92"(zb0.13 
and DEC (J2000)=-66°52'38.2''(±0.9'' 



Manchester et al, 



(2005) then included new and more detailed information 
on the pulsar in a catalogue along with 1532 other pulsars 
and in 2006 re-observed the pulsar at Parkes at a wave- 
length of 20 cm. jFilipovic et al. ( 1998| detected the radio 
source LMC B0530-6655 and suggested it to be an SNR can- 
didate based on its non-thermal spectral index. We note 
that the SNR J0529-6653 (or LMC B0530-6655) lies just on 
the northeast edge of the arc of Ha nebulosity DEM L214 
dDavies et "al] ri976). Finally, [Haberl k Pietsch[ ( [1999 



hereafter HP99) detected a nearby ROSAT X-ray source 
([HP99] 440) at a position of RA (J2000)=05^29^59.3" and 
DEC (J2000)=-66°53'13''. 

Here, we report on radio-continuum. X-ray and optical 
observations of the candidate LMC SNR J0529-6653 with 
its possible association with the LMC pulsar PSRB0529-66. 
Observations, data reduction and imaging techniques are de- 
scribed in Sect. [2] The astrophysical interpretation of the 
moderate-resolution total intensity and polarimetric images 
are discussed in Sect.jS] 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Radio-continuum 

We used radio observations at four frequencies (Table [T]) to 
study and measure flux densities of SNR J0529-6653. For the 
36 cm (Molonglo Synthesis Telescope - MOST) flux density 
measurement given in Table [l] we u sed unpublished images 
as described by [Mills et al. (1984) and for the 20 cm we 
used image from Hughes et al. (2007). Two Australia Tele- 
scope Compact Array (ATCA) projects (C634 and C797; at 
6/3cm) observations were combined with mosaic observa- 



tions from project C918 (Dickel et al. 2005 ). Data for project 



C634 were taken by the ATCA on 1997 August 2, using the 
array configuration EW375. Four days of observations were 
taken from project C797: 1999 May 1-2 (array configura- 
tion: 1.5C), 1999 July 21 (array configuration 750D), and 
1999 July 31 (array configuration: 1.5D). For the final image 
(stokes parameter /) we exclude baselines created with the 
gth p^Qp^ antenna, leaving the remaining five antennas to 
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Figure 1. ATCA image of SNR J0529-6653 at 3 cm overlaid with 
6 cm contours. The contours are from 3 until 17 a with spacings of 
1 cr (0.11 mJy). The black ellipse in the lower left corner represents 
the synthesised beam width of 17.3 x 14.0^^ The circle in the 
centre marks the midpoint of the SNR. The circle in the east 
shows the catalogued X-ray location of [HP99] 440 and the circle 
in the north is the position of the pulsar PSRB0529-66. 



Table 1. Integrated flux densities of SNR J0529-6653. The flux 
density at A=36 cm was estimated using images from [Mills et al.| 
( |1984| > and at A=20 cm from [Hughes et al.| ( |2007t . 
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be arranged in a compact configuration. C634 observations 
were carried out in "snap-shot" mode, totalling ^-^1 hour of 
integration over a 12 hour period. Source PKS B 1934-638 
was used as the primary calibrator and source PKS B0530- 
727 was used as the secondary calibrator. The MIRIAD fSault| 



& Killeen 2006 ) and karma ( Gooch 2006 ) software packages 



were used for reduction and analysis. The 6 and 3 cm images 
(Fig.fT) were constructed using miriad multi- frequency syn- 
thesis ( Sault Wieringa|1994 ) . Deconvolution was achieved 
with the CLEAN and RESTOR tasks with primary beam cor- 
rection applied using the linmos task. Similar procedures 
were used for the U and Q stokes parameters. 



2.2 Optical 

We have examined images (Fig. [2| from the Magellanic 
Clouds Emission Line Survey (MCELS) ( [Smith et al.|2006 ) 
and higher-resolution Ha images (Fig. [s]) obtained with the 
MOSAIC II camera on the Blanco 4-m telescope at the Cerro 
Tololo Inter- American Observatory. The extended arc of the 
H II region DEM L214 is readily seen but we do not find any 
optical nebulosity associated with the SNR candidate or the 
PSRB0529-66. Additional long-slit observations were per- 
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formed on 2010 September 20, using the 1.9-m telescope and 
Cassegrain spectrograph at the South African Astronomi- 
cal Observatory (SAAO) in Sutherland. While the observing 
conditions were good, we did not detect any emission from 
the SNR candidate. 



2.3 X-rays 

ROSAT performed a raster of pointed observations with 
the Position Sensitive Proportional Counter (PSPC) to map 
the soft X-ray emission from the hot gas in the LMC su- 
per giant shell-4 ( Me aburn|[T98Ql [Bomans Dennerl||1999| 
SGS LMC-4). A number of these observations with typical 
exposures of 1000 s include SNR J0529-6653 (Figs, [s] and |4|. 
The HP99 catalogue includes a very weak detection within 
the extent of the radio emission from SNR J0529-6653 which 
does not allow to derive much about the X-ray proper- 
ties of the source (the position of [HP99] 440 is shown in 
Fig. [T]). The hardness ratios are undefined and therefore 
give no information on the spectrum. The catalogue of 
HP99 was derived from the individual PSPC observations, 
not utilising the combined exposure of overlapping fields. 
To investigate [HP99] 440 and its possible association with 
SNRJ0529-6653 and PSRB0529-66 in more detail we se- 
lected 13 observations from the SGS LMC-4 raster which 
covered the SNR/PSR within 24^ of the optical axis (to avoid 
the degraded point spread function at larger off-axis angles) . 
We produced images in different energy bands (broad: 0.1- 

2.4 keV, soft: 0.1-0.4 keV, hardl: 0.5-0.9 keV and hard2: 
0.9-2.0 keV) from the merged data. A colour image of the 
area around the SNR with net exposure (vignetting cor- 
rected) of 11.0 ks is shown in Fig. [4] with red, green and 
blue representing the X-ray intensities in the soft, hardl 
and hard2 bands. The resolution of the ROSAT PSPC varies 
with energy but the point spread function is always less than 

v. 



3 RESULTS AND DISCUSSION 

The radio-continuum remnant SNR J0529-6653 exhibits 
a ring-like morphology, indicative of a shell struc- 
ture, with brightened region along the southwest rim 
(Fig. [l|. It is centred at RA (J2000)=05^29^51.0" and 
DEC (J2000)=-66°53'27.l''. We estimate the spatial extent 
of SNRJ0529-6653 (Fig. ^ at the 3cr (Table [l] Col. 4) 
level (0.33 mJy) along the major (N-S) and minor (E-W) 
axes. Its size at 6 cm is 137'' x 128''=b4'' (33x31 pc with 
1 pc uncertainty in each direction). We also estimate the 
SNRJ0529-6653 ring thickness to <20'' (5 pc) at 6 cm, 
about 30% of the SNR's radius. 

An extended X-ray source is clearly seen at the loca- 
tion of the SNR with a brighter spot right at the position 
of [HP99] 440. Our images are based on the combination of 
data from 13 observations, which had too short exposures 
to allow an astrometric alignment to a common reference 
system. Therefore, in order to verify the source extent we 
compared the projected profiles of the extended source and 
other nearby sources, which appear point-like. For the pro- 
file, we projected the counts in the 0.1—2.4 keV image within 
a slice of 220'' width and position angle of 330° (to investi- 
gate the profile of the bright spot). We measure a FWHM of 





[Sii] . 




RA (J2000) 



30^ 



Figure 2. MCELS Ho; {top), [Sll] {middle) and [Olll] {bottom) 
images of the area around SNR J0529-6653 overlaid with 6 cm 3a 
(0.33 mJy) radio-continuum contour. 
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Figure 3. The MCELS-2 Ha image of the area SNR J0529-6653 
overlaid with the ROSAT PSPC contours from the broad-band 
(0.1-2.4 keV) image (pixel size: 15^^). Contours are from 3 to 
7 cts/pixel with 0.5 cts/pixel spacings. 



130^^=b7^^ for the extended source and 75^^±7^^ for the source 
visible at the southern rim of the image in Fig. |4] which ex- 
hibits a similar peak intensity. The extended profile shows a 
peak on top of the broader intensity distribution confirming 
the brightest part of the X-ray emission in the north-east. 
The presence of extended X-ray emission is coincident with 
the radio emission of the SNR (Fig. [g]). 

The latter is also supported from the relation between 
X-ray luminosity and energ y loss for rotation- powere d 
pulsars (L 



=(10- 

With a value of E 



MO-^)xE; 
=6.6x10^2 



Becker & Truemper (1997)). 

T 



erg s " (from online ATNF pul- 
sar catalogu^ for PSR B0529-66 we expect an X-ray lumi- 
nosity for the pulsar of well below 10^^ erg s~^ , several orders 
of magnitude below the detection limit of ^ 10^^ erg s~^ for 
a 10 ks ROSAT PSPC observation. Also, from the estimated 
spin down age of ^10^ yr (see also below) one would not ex- 
pect the existence of a PWN. 

The net counts of 337=b25 cts s~^ derived from the SNR 
in the 0.1—2.4 keV band are insufficient for spectral mod- 
elling in order to derive temperature and density from which 
one could estimate the explosion energy. Assuming a thermal 
emission model (APEC) with a temperature kT=0.19 keV, 



typical for faint SNR in the Magellanic Clouds (Filipovic 
et al.|2008 ) and a foreground absorption of 6x10^ H cm~ 
we can roughly estimate the absorption-corrected X-ray lu- 
minosity to 5.7x10^^ erg s~-^. 

We estimate the radio spectral index a=— 0.68±0.03 
from the integrated flux densities given in Table 1. This 
spectral index is consistent with those of known SNRs of 
mid-to- younger ages. SNR J0529-6653 is located in the cen- 
tre of the SGS LMC-4 (some -3° from 30 Dor) and there- 
fore, conceivable that it is expanding in a very low density 



http: / / www.atnf.csiro.au / research / pulsar / psrcat / 
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Figure 4. The ROSAT PSPC RGB colour image of the area 
around SNR J0529-6653. The energy bands are: red (0.1-0.4 
keV), green (0.5-0.9 keV) and blue (0.9-2.0 keV). The image 
has a pixel size of 15^' and is smoothed with a cr of one pixel. 
A white ellipse is centred on the position of the SNR J0529-6653 
with extent of 137^^x128^^ (33x31 pc). 




Right Ascension (J2000) 

Figure 6. ATCA image of SNR J0529-6653 at 6 cm overlaid with 
the contours from ROSAT PSPC image (0.1-2.4 keV). Contours 
are as in Fig. [3] The ellipse in the bottom left corner represent 
the beam of 17.3'' x 14.0'^ 



environment, which causes a slightly steeper spectral index 
for its age (see Sect.js] Para. 10). 

Fig. [t] shows a surface brightness-diameter — D) 
diagram at 1 GHz with theoretically-derived evolutionary 
tracks ( |Berezhko k Volk|2004t superposed. SNR J0529-6653 
lies at (D,!]) = (32 pc, 2.3 x 10"^^ W m-^ Hz"^ Sr"^) on 
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Figure 5. 

I-Profile of 



The left image shows the 6 cm image overlaid with the major and minor axis labels. The centre and right images show the 
the major and minor axis respectively, with an overlaid black line at 3cr. 



the diagram. Its position tentatively suggests that it is in 
the early Sedov phase of evolution — expanding into a very 
low density environment with the higher initial energy of 
a supernova explosion (2 — 3 x 10^^ ergs) and the age of 
^25 000 yr. We acknowledge that such scenario couldn't ex- 
plain for the lack of SNR's optical emission. Alternatively, 
this SNR could be bit older (^70 000 yr; the end of adiabatic 
phase) and expanding into the medium dense environment 
with the minimal start-up energy. 

A linear polarisation image was created for 6 cm wave- 
length using Q and U parameters (Fig. [8|. The mean frac- 
tional polarisation was calculated using flux density and po- 

\/ ^'q'^ ^'u 

larisation: P=-^^^ — • 100%, where Sq,Su and Si are 

integrated intensities for the Q, U and / Stokes parameters. 
Our estimated peak value is 17%±7% {3a) at 6 cm and no 
reliable detection at 3 cm. Along the south side of the SNR 
shell there is a pocket of uniform polarisation (Fig.[8|, possi- 
bly indicating varied dynamics along the shell. The polarisa- 
tion appear to be uniformly distributed (tangential orienta- 
tion) coinciding with the SNR total intensity. As the whole 
SNR is nearly round the only irregularity is the brightness 
variations. While the distance from the nearby Ka filament 
(DEML214) seems bit far away it is most likely that the 
SNR interacted with some less obvious local clouds so the 
radio-continuum shape has not yet been disturbed. 




0.1 1 10 100 



Without reliable polarisation measurements at a sec- 
ond frequency (3 cm) we cannot determine the Faraday ro- 
tation and thus cannot deduce the magnetic field strength. 
However, by using the new equipartition formula for SNRs 
( Arbutina et al.| |2011), we can estimate the magnetic field 
strength for the SNR J0529-6653. The derivation of the new 



equipartition formula is based on the Bell (1978) diffuse 



shock acceleration (DSA) theory. This derivation is purely 
analytical, accommodated especially for the estimation of 
magnetic field strength in SNRs. Using this new formula, 
the calculated magnetic field strengths for SNRs are be- 



tween those calculated by using classical equipartition (Pa- 



cholczyk 1970) and revised equipartition (Beck & Krause 



2005 ) . The average equipartition field over the whole shell of 



D, pc 



Figure 7. Surface brightness-to-diameter diagram from |Berezhko| 
|fc Volkl ( [2004t , with SNR J0529-6653 added. The evolutionary 
tracks are for ISM densities of ^h= 3, 0.3 and 0.003 cm"*^ and 
explosion energies of F^sn = 0.25, 1 and 3x10^-'^ erg. 

corresponding "calculator"]^, corresponding those of young 
to middle-aged SNRs where the interstellar magnetic field is 
compressed and amplified by the strong shocks. 

The Parkes pulsar catalogue suggests that the pul- 
sar PSRB0529-66 is 9.97± 0.77x10^ yr old. SNRs nor- 



The 



calculator 



ilable 



SNR J0529-6653 is ^48 /iG (see Arbutina et al. (2011); and http://poincare.matf.bg.ac. rs/~arbo/eqp/| 
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mally merge into the ISM and become unrecognisable af- 
ter ^150 000 yr. We point out that a measured radius of 
^-^17 pc would imply a mean expansion speed of only about 
0.1 km s~^ for a million year age which is clearly unreason- 
able. Therefore, this factor of ^6 in age difference would rule 
out the possibility that the pulsar and SNR were formed 
from a common supernova explosion. However, IMigliazzo] 
et al. ( 2002 ) showed that pulsar age estimates from its "char- 



limit. However, unless the direction of the pulsar motion is 
significantly out of the sky plane, its real velocity will not 



acteristic age" (tc = P/2P) could be uncertain when com- 
pared to its corresponding SNR age. Also for nearby Galactic 
isolated neutron stars large differences between the charac- 
teristic age and the dynamic age estimated from their proper 
motion and the likely birth place are found (a factor of ~8 
for the case of RX J1856. 5-3754, see [Tetzlaff et al.|[20TT ). 
They also argue that the ratio depends on the magnetic 
field strength and its evolution in the past. 

If we assume that the SNR J0529-6653 and pulsar 
PSRB0529-66 are associated and that the pulsar was born 
near the geometric centre of the SNR, then the pulsar 
must have moved ~14zbl pc to its current location (Craw- 



ford et al. 2001). Assuming the above "characteristic pul- 



sar age", the pulsar would be travelling (in the sky plane) 
at ~ 14 km s~^, an absolute minimum kick velocity. If we 
assume a canonical SNR age of 25 000-50 000 yr then the 
pulsar velocity in the sky plane would be in a range of 550 
to 275 km s~^, in agreement with the typical pulsar kick ve- 
locities (| Lyne Lorimer|1994 ). The direction of the pulsar's 
travel is not necessarily in the sky plane, so this is a lower 



be much greater. Klinger et al. (2002) found a jet-like struc- 
ture and a point X-ray source within t he LMC SNR N 206 
travels at a speed of 800 km and Owen et al. (2011) 



estimated the kick velocity of the pulsar candidate inside 
SMC PWN IKT16 to be approximately 580 km s"^ The 
possible pulsars have left prominent radio trails in N 206 and 
IKT 16, but no such trail can be seen in our radio images of 
SNR J0529-6653. This makes a SNR-pulsar connection less 
likely. Nevertheless, if PSRB0529-66 is related to the SNR, 
it is likely to be travelling more slowly and thus any trail 
would be less prominent. We acknowledge that even if the 
pulsar is not traveling, there may still be a prominent PWN 
in the X-rays, although it will not be trailed. Unfortunately, 
the lack of evidence of PWN in the ROSAT data is purely 
because of resolution, and therefore, cannot be taken as an 
evidence of lacking a PWN. 

The absence of detectable optical emission (Fig.[2| from 
this SNR is not unique as other well-studied SNRs, such as 
LMC SNR J 0528-6714 ([Crawford et al.|2010| or the Galac- 
tic Vela Jr. (Stupar et al. 2005), also do not exhibit opti- 
cal emission. The ISM in the interior of SGS LMC-4 must 
have very low density and hence the lack of optical emission. 
We argue that, as with the majority of other SNRs in the 
Magellanic Clouds, this intriguing SNR is most likely in the 
adiabatic phase of its evolution ( Payne et al.|2008J simply 
because of its modest size. 
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In the MCELS-2 Ha image (Fig. [3| we clearly see the 
large extent of Hii region DEML214 in the south, south- 
west and south-east. We also see smaller scale filamentary 
structures in the north of the remnant which, together with 
DEM L214 and the Ha emission around the cluster in the 
east, seem to form a large scale shell. This and the distribu- 
tion of the massive stars (see next paragraph) indicate that 
some heating must have already taken place where the SNR 
is located. A shock expanding in a already hot thin medium 
is not efficient and expands quickly. This may explain the 
faint X-ray appearance of the remnant SNR J0529-6653. 

As massive stars rarely form in isolated environments, 
core-collapse supernovae are most likely superposed on a 
stellar population rich in massive stars. We have used 
the Magellanic Cloud Photometric Survey ( [Zaritsky et al. 
|2004( MCPS) data to construct colour-magnitude diagrams 
(CMDs) and identify blue stars more massive than ^8 Mq 
within the area of SNR J0529-6653. The CMD in Fig. [5] 
(left) indicates that most of the blue massive stars are main 
sequence B stars. Their spatial distribution marked in blue 
Fig. [9] (right) shows higher concentrations toward a cluster 
to the east and the nebulosity to the southeast of the SNR, 
but not within the SNR and its immediate vicinity. The pro- 
genitor of SNR J0529-6653 could be a B star (core-collapse 
supernova) or an accreting white dwarf in a binary system 
(Type la supernova) . To distinguish between these two pos- 
sibilities, deep XMM-Newton X-ray observations with suffi- 
cient counts for spectral analysis of plasma abundances are 
needed. 



4 CONCLUSION 

We have carried out the first detailed multi-frequency study 
on a recently-detected LMC SNR J0529-6653, which pre- 
viously had records for a pulsar at this position. We esti- 
mated a diameter of 137'' x 128''=b4'' (33x31 pc with 1 pc 
uncertainty in each direction), a spectral index {a = 
-0.68±0.03) and a peak polarisation of -17%±7%. We note 
that even though there is no optical emission (and thus no 
[S ll]/Ha>0.4), the presence of the non-thermal radio and 
the X-rays do satisfy two of the three criteria for calling 
it an SNR. While there is a positional association between 
the SNR and PSR, they are not necessarily related as the 
ages are inconsistent. While assuming that the PSR has the 
same age as the SNR yields estimated kick velocities that 
are fairly typical, there is no radio-frequency trail of the kind 
seen in other SNR-PSR connections. 
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Figure 9. Left panel shows a B — y vs. V colour-magnitude diagram from the MCPS ( [Zaritsky et al.|2004| ). Stellar evolutionary tracks 
from |Lejeune S chaerer (2001) are shown with solid red lines (5, 15, 25 and 60 Mq) and dashed red lines (3, 9, 20 and 40 Mq). The 
selection criteria used for probable B-star candidates is shown as a heavy dashed blue lines. On the right is the MCELS-2 Ho; image of 
the area around SNR J0529-6653 overlaid with probable B-star candidates (blue crosses) from the MCPS. A 90^^ radius circle (red) is 
centred on the position of the SNR. 
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